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ABSTRACT 



Aims. We report on the first spectroscopic search for quiescent degenerate binaries in globular clusters. 

Methods. Our survey is based on a sample of short-period optical variables which are likely optical counterparts of quiescent X-ray 
sources in cj Centauri (NGC 5139) and NGC 6397. 

Results. The studied candidates have nearly sinusoidal light curves with amplitudes of 0.05-0.12 mag in V (0.35 mag in one case), 
and periods of 0. 1-1 .3 days. This type of variability, most probably originating from the ellipsoidal effect, has been observed in X-ray 
novae when they settled into quiescence after an outburst. 

Conclusions. We find that two of the surveyed systems harbour dim components with masses in excess of 1 M G , making them 
attractive targets for future investigations. We also suggest that there are two subpopulations of blue stragglers in lj Centauri, differing 
in mass-transfer history and/or helium content. 
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stragglers 
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1. Introduction 

Since the early surveys with the Uhuru and OSO-7 satellites 
it has been known that X-ray sources occur about a thousand 
times more frequent ly in globular clusters (G Cs) than in the rest 
of the Galaxy (e.g., lKatzlll975b IClarklll975l) . Over fifteen hun- 
dred of them have been detected to date in these stellar sys- 
tems by recen t missions such as Chandra and XMM-Newton 
(Poolev 2010, and references therein). In addition, over a hun- 
dred millisecond pulsars - objects which have evo lved in low- 
mass X-ray binaries - are known to reside in GCs dLvnch et al.l 
120101) . Most of the X-ray sources are binaries with degenerate 
components (neutron stars or white dwarfs), henceforth referred 
to as "degenerate binaries". Remarkably, no stellar-mass black 
holes (BH) are known in Galactic G Cs, despite theoretica l pre- 
dictions that they should exist (e.g., iDevecchi et al.ll2007h . The 
presence of intermediate-mass black holes (IMBHs) (100M o < 
M < 10 3 ~ 4 tf n ) in G Cs h as also not yet been proven (e.g ., 
M accarone & Servillatll2008l l van der Marel & Andersonl2010l) . 



While the origin of IMBHs is unclear, stellar-mass black 
holes should form along with neutron stars via supernova explo- 
sions during the early evolution of a cluster. The very presence of 
such objects would provide interesting information concerning 
the effectiveness of dynamical interactions leading to the ejec- 
tion of BHs into the intracluster medium, which according to 
the predic tions of current appro ximate models should be very 
high (e.g.. iDowning et a"fll2010l) . Their absence would confirm 
those predictions. Alternatively, it could mean that massive low- 
metallicity stars cannot produce BHs, which would have impor- 
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tant implications for the origin of early BHs, believed to be seeds 
for the first galaxies in the Universe. In either case, searching for 
black holes in GCs is certainly a worthwhile task. 

It is well known that field X-ray novae (binaries most proba- 
bly hosting stellar-mass black holes) spend most of their time 
in quiescenc e, showing only ellipsoidal vari ability in the vis- 
ible domain dRemillard & McClmtockll2006l) . Since GCs har- 
bor a multitude of degenerate binaries, they must also contain 
many systems which presently accrete at a very low rate or 
do not accrete at all, thus being very weak or undetectable in 
the X-rays. Based on X-ray luminosity (Lx) and hardness ra- 
tio, Grindlay (2006) splits the population of weak X-ray sources 
in GCs into four major classes. The first three, arranged ac- 
cording to Lx increasing from 10 29 (the present sensitivity limit 
for the nearest GCs) to 10 32 erg s _1 are i) active binaries (bi- 
nary main-sequence stars in which the X-ray emission originates 
from chromospheric activity, e.g. BY Dra type systems); ii) cat- 
aclysmic variables (in which a white dwarf accretes from a low- 
mass main-sequence companion); iii) quiescent low-mass X-ray 
binaries (in which a neutron star intermittently accretes from a 
main-sequence or evolved companion), hereafter referred to as 
qLMXBs. The fourth class is composed of millisecond pulsars, 
whose luminosity strongly depends on the predominant mech- 
anism generating X-ray quanta (residual accretion; collision of 
the pulsar wind with the ambient medium or with matter lost by 
the companion; thermal surface emission). 

The study of weak X-ray sources is hampered by the fact that 
even in the nearest GCs they are detected down to the present 
sensitivity limit, and often remain unclassified as they are too 
dim to allow for an estimate of the hardness ratio. Most of them 
are likely active binaries or cataclysmic variables, but some may 
be quiescent systems with neutron stars or even black holes. 



2 



Rozyczka et al.: A search for degenerate binaries in u> Centauri and NGC 6397 



Optical counterparts of these systems are often weak and hard 
to ide ntify in the crowded GC environment (e.g.. IVerbunt et all 
l2008h . So far, the X-ray data have allowed us to identify a few 
candidate qLMXBs, but their nature has not been confirmed by 
optical spectroscopy dGuillot et al. 2009[). 

We decided to follow an entirely different approach. Instead 
of looking for optical counterparts of X-ray sources, we select a 
sample of candidate systems for degenerate binaries from short- 
period, low-amplitude optical variables catalogued in existing 
surveys, and measure their radial velocities. Our targets have 
nearly sinusoidal light curves with V-band amplitudes smaller 
than 0.35 mag and periods shorter than ~1.3 days. This type of 
variability is common in all classes of close degenerate bina- 
ries; in particular, it is observed in optical counterparts of sev- 
eral Galactic X-ray novae whic h almost certainly harbor BHs 
dRemillard & McCIintockll2006l) . The optical modulation is in- 
duced mostly by the ellipsoidal effect from the nondegenerate 
component, with an amplitude depending mainly on the Roche 
lobe filling factor and the inclination of the orbit. Close degener- 
ate binaries can be unambiguously identified as single-line sys- 
tems with large orbital velocities (K > 150 km s~'). Among 
them, those harboring a black hole are distinguished by a mass 
function /,„ = (m sin 3 /( m 6m) 2 > 2M , where tnun an d m stand 
for the t otal mass of the binary and the m ass of one of its com- 
ponents dRemillard & McClintock1l2006l) . 

Our sample and observational data are introduced in Sect. 2. 
A detailed analysis of the data is reported in Sect. 3, and the 
results are discussed in Sect. 4. 
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Fig. 1. Color-magnitu de diagram of u> Cen based on CASE data 
dKaluznv et alj 120041) trans formed t o the BV system with the 
help of standards listed by IStetsonl d2000l) . Marked are the lo- 
cations of the investigated systems. 
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2. Observations and data reduction 

Optical counterparts of active X-ray sources in GCs have been 
found essentially everywhere in color-magnitude diagrams: to 
the left of the main sequence, on the main sequence, to the right 
of it, on subgiant, giant and horizontal branches, and also above 
the main-sequence turnoff in the regions occupied by blue strag- 
glers and EHB stars (e.g. Heinke et al. 2005; Servillat et al. 
2008). Thus, while preparing this study we felt it justified to 
pick photometrically suitable targets without paying attention to 
their location with respect to the main sequence. Since large- 
amplitude variables are more likely to be ordinary contact bina- 
ries, we avoided objects with full V-band amplitudes in excess 
of ~0.35 mag. Suspected pulsating variables were excluded; we 
also took care to exclude systems with appreciable X-ray emis- 
sion. Where it was possible, the memb ership status of the tar- 
get was established based on the data o f lBeUini et aLld2009h for 
ijj Cen, and Anderson (priv. comm.) and IStrickler et al 1 120091) 
for NGC 6397. In the remaining cases membership was decided 
a posteriori based on the systemic velocity with respect to the 
cluster. 

Alltogether, we selected seven objects in u> Cen and four in 
NGC 6397. In color-magnitude diagrams two of them are lo- 
cated to the left of the main sequence, seven above the main- 
sequence turnoff, and two at the turnoff itself (see Figs.Q]and[2]i. 
A summary of their basic data i s given in Tabled] (the labels are 
the same as those introdu ced bv lKaluznv et al.l d2004l) for a> Cen 
and lKaluznvetal.1 d2006l) for NGC 6397). 

Our pape r is ba s ed on photometric measurements of 
120061) . supplemented by those of 
and by our unpublished data for NGC 
du Pont telescope at Las Campanas 



Kaluznvetal. 
Weldrake et all d20O 
6397 obtained on the 
Observatory, Chile (a total of 771 V-band frames was acquired 
betw een 2009 June 20 and June 30, and reduced in the same way 
as in iKaluznv et al.l d2006l) ). The newer data demonstrated that 



the observed light variations were coherent, although in some 
cases the light curve was changing from symmetric to asymmet- 
ric or vice versa. 

Spectroscopic data were collected during the nights 20/21 
and 21/22 of May 2009 with the MagE (Magellan Echellette) 
spectrograph attached to the 6.5-m Magellan-Clay telescope at 
Las Campanas Observatory. The seeing varied betweeen 075 and 
078 on the first night, and betweeen 074 and 079 on the second 
one. A 0785 slit was used, providing a resolution R = 4820. 

During the observations pairs of scientific spectra taken for 
the same target were separated by an exposure of a thorium- 
argon hollow-cathode lamp. The exposure times per spectrum 
ranged from 120 s to 1080 s, depending on object brightness and 
observing conditions. After bias and flat-field correction each 
pair of frames was combined into a single frame, allowing for 
the rejection of cosmic ray hits. The observations were reduced 
with the IRAlQ ECHELLE package. 

In order to determine the mass function of a given target it 
was sufficient to take just a few spectra at different phases, and fit 
a simple sinusoid to its velocity curve phased with the photomet- 
ric ephemeris. Altogether, 96 spectra were obtained (from 5 to 
10 per object). The useful range of the reduced spectra extended 
from 4050 A to 6800 A, where the spectra had 13 < S/N < 23 
(with a few cases of lower quality). Radial velocities were mea- 
sured in that range with the help of IRAF routines FXCOR and 
XCSAO, using synt hetic templates from the library compiled by 
iMunari et al ] d2005l) . The results were verif ied by applying the 
broadening function formalism described bv lRucinskl d2002l) . In 
most cases the differences between various measurements were 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the AURA, Inc., under co- 
operative agreement with the NSF. 
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Cluster Star a 200 o 62000 Pm* L x ** V B-V ref.f 
[h:m:s] [":':"] [10 30 erg s"'] [mag] [mag] 



a> Cen 


V214 


13:27:21.82 


-47:37:19.0 


97 


<13 


16.55 


0.25 


1,2 


a> Cen 


V240 


13:27:28.68 


-47:26:19.5 


99 


<L2 


16.31 


0.30 


1.2 


a> Cen 


V251 


13:27:28.02 


-47:26:43.7 


97 


<1.2 


16.17 


0.25 


1 


oj Cen 


V254 


13:27:28.61 


-47:27:39.0 


99 


<1.2 


16.41 


0.25 


1,2 


a> Cen 


NV334 


13:27:23.39 


-47:29:11.9 


88 


<1.2 


18.63 


0.35 


1 


a> Cen 


NV360 


13:26:02.06 


-47:32:23.6 


100 


<1.2 


17.09 


0.25 


1 


a> Cen 


NV400 


13:26:10.09 


-47:31:50.3 


99 


<1.2 


16.56 


0.20 


1 


NGC 6397 


V17 


17:40:43.78 


-53:41:16.2 


100 


0.4 


16.17 


0.40 


3,4 


NGC 6397 


V20 


17:40:41.66 


-53:40:33.1 


100 


1.6 


15.75 


0.40 


3.4 


NGC 6397 


V27 


17:41:13.80 


-53:41:14.1 






18.19 


0.45 


3.4 


NGC 6397 


V36 


17:40:44.10 


-53:42:11.3 




0.4 


16.48 


0.60 


3,4 



* Membership probability: [Bellini et all d2009h for V214-NV400: IStrickler etal] d2009h and Anderson (2010, priv. comm.) for V17-V20; no data 
available for V27 and V36. 

** X-ray luminosity: iGendre et all d2003h for V214; lHaggard et all d2004h for V240-NV400; iBogdanov etal] d2010|) for V17-V20 and V36; no 
data available for V27. 

t Photometry: 1 - iKaluznv etal] d2004h : 2 - IWeldrake eTal] d2007l) : 3 - iKaluznv etai] d2006l) : 4 - our unpublished data collected in 2009. 
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Fig. 2. Color-magn itude diagram of NGC 6397 (after 
IKaluznv et al.l 120061) with marked positions of the investi- 
gated systems. 



smaller than Av = 10 km s , which may be regarded as a fair 
estimate of the error. 



3. Analysis 

While our data are insufficient for a detailed analysis, an esti- 
mate of the basic parameters of the observed systems is possi- 
ble. In all cases but one (V20 in NGC6307) single-line spectra 
are observed, indicating that to a first approximation all the light 
we receive from the binary is produced by just one of the com- 
ponents. Based on this assumption, we may set some limits for 
acceptable configurations. 



3.1. Method 

Our observables are: orbital period P, apparent V-band magni- 
tude V, color index B-V, amplitude of the light curve AV bs> an d 
amplitude of the velocity curve K b s . Since the distance moduli 
to a> Cen and NGC 6397 are known, V of a system belonging 
to the cluster is equivalent to the observed bolometric magni- 
tude of its visible component M£ o/ ("observed" means here "de- 
rived from observations; hence index "o"). We a dopt (m - M)p = 
13.75 mag and E(B - V) — 0.1 15 mag for u> Cen dVillanova et al.l 
l2007l and references therein), and (m — M)g - 12.03 mag and 
E(B - V) = 0.183 mag for NGC 6397 (iRicher et alj|2008f and 
references therein). Assuming Ay = 3.1 E(B - V), the corre- 
sponding values of Ay are 0.36 mag and 0.57 mag, yielding a dis- 
tance modulus (m-M)v of 14.1 1 mag for ai Cen and 12.60 mag 
for NGC 6397. 

The effective temperature of the visible component T is 
known from the relation between the color index B-V and 
temperature derived for an appropriat e chemical compos ition. 
We adopt [Fe/H] = -1.8 for oj Cen (ISollima et al] |2009l) and 
[Fe/H] = -2.0 for NGC 6397 (ICarretta et al]|2009l). and we us e 
color transformation tables provided by Leieune et al.l dl998l) . 
The parameters derived for our target stars from the observa- 
tional data are given in Table [2] 

The last column in Table [2] contains mass-center velocities 
with respect to the cluster obtained by a simple sinusoid fit, 
with heliocentric velocities of the clusters (232.3 km s~' for 
lo Cen and 18.9 k m s' 1 for NG C 6397) taken from the cata- 
logue compiled by lHarrisI i 19961) and updated online The ve- 
locities of V214, V240, and V254 may seem a bit largish for 
their distances from the cluster center (17'0, 14'9, and 14'7, 
respectivel y), even though the published cluster tidal radius is 
r, = 44f8 (irageret all ll 9951) . or r, = 57:03 according to the 



1 ~ TT v ^ a ^ L ^ aL " 177 -"'- U1 't ~ a^uiumg iu Lilt 

lHarrisI i 19961) catalog. That these stars are well wi thin the limits 
for clu ster membership is confirmed by Fig. 6 in ISollima et al.l 
(120091) . which shows the variation in radial velocity of u> Cen 
stars as a function of distance from the cluster center. Still, V240 
turned out to be a doubtful case (see Sect. 13.3. 11 1. As for V214 
and V254, we checked that the best fits described in Sect. 13.3.11 
were not affected when vq was lowered, respectively, to 10 and 



http://www.physics.mcmaster.ca/Globular.html 
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Table 2. Parameters of the target objects derived from the observations. 



Cluster 


Star 


P 


M bol 


T 


K-obs 


fm 


v 






[d] 


[mag] 


[K] 


[km s 1 ] 




n -In 

[km s '] 


a> Cen 


V214 


0.341806 


2.44 


8300 


35 


1.5xl0~ 3 


21 


co Cen 


V240 


0.331888 


2.19 


8000 


< 10 


< 3.4 x 10~ 5 


21 


a> Cen 


V251 


0.922458 


2.06 


8300 


35 


4.1xl0~ 3 


6 


co Cen 


V254 


0.385071 


2.30 


8300 


60 


8.0x10 


18 


co Cen 


NV334 


0.257881 


4.52 


7300 


35 


l.lxlO" 3 


-9 


co Cen 


NV360 


0.630821 


2.98 


8300 


45 


5.9xl0~ 3 


-17 


co Cen 


NV400 


0.636874 


2.40 


8800 


20 


5.3xl0~ 4 


-7 


NGC 6397 


V17 


1.061316 


3.35 


6200 


50 


1.4xl0- 2 


-2 


NGC 6397 


V20 


0.861177 


2.94 


7200 


* 


* 


9 


NGC 6397 


V27 


0.556134 


5.38 


7500 


25 


9.0xl0" 4 


-6 


NGC 6397 


V36 


1.098569 


3.65 


6200 


95 


9.7xl0~ 2 


-1 



* Velocities of both components were measured, yielding K\ = 25 km s 1 and Ki = 130 km s . The corresponding mass functions are f„, 
1.4 x 10~ 3 M Q and f m2 = 1.9 x 10~' M . 



6 km s -1 . The velocity of NV360, although also rather large, 
is entirely consistent with the distance of this system from the 
center of u> Cen (3 '6), at w hich the relative velo cities of cluster 
members reach +35 km s _1 dSolrima et a l. 2009). 

With T being fixed, the task is to find the orbital separation 
a, inclination of the orbit i, and masses of the components ni\ 
and »■ Unless indicated otherwise, we assume that the system 
is semi-detached, with the visible component filling its Roche 
lobe. It must be stressed here that this particular model of the 
binary is not related to the physical structure of our targets (in 
fact, in some cases it might imply a mass transfer unstable on 
a dynamical scale). A semi-detached system is simply the most 
compact one among those with Mf , = M? ,. Focusing on it, we 
minimize the orbital separation, so that the corresponding "semi- 
detached" masses mf and mf are the smallest allowable given 
P and q. The lower limits of the actual masses of the compo- 
nents are obtained by finding "the best" mass ratio q/„ defined 
as the one for which AV c /AV j, s = 1 (see Sect. 13.21 for further 
discussion). Our analysis involves the following steps: 

1 . Specify if the visible component is primary or secondary. 

2. Set the mass ratio q = milm\ < 1 (we define the primary as 
the more massive component, not necessarily being the more 
luminous one). 

3. Given the orbital period P, adjust the separation a so that the 
calculated bolometric magnitude M c bo[ of the visible compo- 
nent be equal to M° boV 

4. Adjust i and vo so that the calculated velocity curve fits the 
observed one, and find the ratio of the computed amplitude 
of the light curve (AV C ) to the observed amplitude (AV b s ). 

5. Repeat steps (2)-(4) for several q values. 

6. If in step (1) the visible component was specified as the pri- 
mary, specify it as the secondary (and vice-versa). Repeat 
steps (2)-(4) for several values of q. 

The calculations in ste ps (3) and (4) are p erformed using 
the PHOEBE i nterface ilPrsa & Zwitterl 120051) to the Wilson- 
Devinney code (I Wilson & Devinnevl 19711) . 

In principle, u sing a relation betwe en V, B - V and stellar an- 
gular diameter dKervella et al. 1 120041) would be more straight- 
forward and simpler than calculating M°, , and T, and comparing 
M° bol to M c boI (note that T is needed for PHOEBE to derive M € bol ). 
Based on such a relation, we could directly fit the "observed" ra- 
dius of the visible component R° = QD (where D is the distance 
to the cluster) to the calculated average radius of the Roche lobe. 



Unfortunately, fits of iKervella et alj d2004l) proved to be unreli- 
able when extrapolated to low metallicity and small 6. This is 
not surprising, as these authors explicitly warn about the nonlin- 
earity of the relation involving V and B - V. 



3.2. Discussion of errors 

The values of our key input parameters, i.e. temperatures and 
bolometric luminosities, are not known accurately, and their er- 
rors are difficult to estimate without engaging in a lengthly (and 
likely ambiguous) discussion of involved factors. In order to ver- 
ify how these uncertainties influence the output we varied T and 
M° bgl by small amounts and observed the corresponding varia- 
tions in m\ and mo. We found that decreasing T by 100 K causes 
the masses to grow by 7 
-0.1 mag makes them larger by 10 - 15% 

At a first glance, focusing on semi-detached systems seems 
too restrictive. To see it clearly, suppose that the visible compo- 
nent is entirely contained within its Roche lobe. Then in order to 
recover M° bol we would have to increase the orbital separation a. 
Since the period P is fixed, the masses would have to increase, 
too. To keep the observed velocity amplitude constant, the in- 
clination ; would have to decrease. For that reason, and also be- 
cause the visible component would now be less deformed, AV C 
would also decrease. This way the range of q for which the ra- 
tio AV c /AV bs ~ 1 could be extended onto fits which for semi- 
detached systems yield AV c /AV b s > 1. In principle, it is even 
possible that mi and mo derived for some mass ratio q + qb 



10%, while an increase in M° bol by 



from detached fits could be smaller than mf{qb) and mf(qb). 
However, upon verifying this possibility we found that in most 
cases the detached fit results in either too low an amplitude for 
the light-curve or masses larger than those obtained for qb, so 
that the extension of the g-range is marginal (if any). We con- 
clude that the errors in the mass limits caused by neglecting de- 
tached configurations are smaller than those related to bolomet- 
ric corrections or the color- temperature scale. 

Another possible source of errors is the assumption that the 
whole light output of the system originates in the component 
responsible for the observed spectrum. The components of the 
only system in our sample in which two sets of velocities could 
be measured differ in brightness by AM « 1 mag. Assuming 
that in the remaining systems AM = 2 mag causes their brighter 
components to be dimmer by 0.16 mag. The corresponding de- 
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crease in masses, estimated from the general relation between 
mass and M? , mentioned earlier, amounts to 15 - 20%. 

3.3. Results 

Quantitative results of the analysis are shown in Tables 3-13, 
and discussed in Sect. 13.31 The best fits are indicated with an 
asterisk in the last column of each table, and the corresponding 
computed light and velocity curves are displayed in Figs. [3j |4] [5] 
and [6] along with the observational data. 

3.3.1. cj Centauri 

V214. Assuming that the visible component is the secondary re- 
sults in computed light-curve amplitudes consistently much too 
low. Assuming that it is the primary makes the difference be- 
tween AVobs an d AV C smaller, but still unacceptable. 

In order to make AV C larger, we are forced to look for solu- 
tions in which both components contribute to the observed ellip- 
ticity effect. Based on the simplest possible assumption that the 
system is in contact, we apply the following procedure: i) set q\ 

ii) adjust a so that the combined M\, i of the system is 2.44 mag; 

iii) adjust i so that the computed velocity curve agrees with the 
observed one; iv) check if the computed light curve agrees with 
the observed one. 



Table 3. Derived parameters for V214 



9 


a 


i 


mj 


mi 


AVc/AVot, 




[Re] 


n 


[M ] 


[M G ] 




0.10 


2.45 


62 


1.54 


0.15 


1.1 


0.13 


2.55 


58 


1.69 


0.22 


1.0* 


0.15 


2.60 


47 


1.76 


0.26 


0.6 



The fits are shown in Table [3] and we conclude that the 
available data favor a system with a rather large mass ratio, 
which must have undergone significant mass transfer (and 
may still be transferring mass at a low rate). Obviously, the 
sum m = mi + '«2 should not exceed two t urnoff masse s of the 
youngest population of at Cen, i.e. 1.84 M (lNorrisl2004l) . Given 
the rather large errors in the mass estimates due to uncertainties 
in T and Mf ,, one may accept that this requirement is fulfilled 
by at least the first two fits in Table [3] 

V251. Like in the previous cases, assuming that the spectrum 
originates in the secondary we find that for all q values AV C 
is much smaller than AV b S . Assigning the spectrum to the 
primary leads to the results displayed in Table [4] As the total 
mass is markedly smaller than 1.84 M Q , we seem to have a 
system which not only underwent mass transfer, but also lost a 
significant amount of mass. 



Table 4. Derived parameters for V25 1 



q 


a 


i 


m l 


m.2 


AVJAV obs 




[Re] 


[°] 


[M e ] 


We] 




0.20 


3.10 


33 


0.39 


0.08 


3.2 


0.35 


3.44 


32 


0.48 


0.17 


1.0* 


0.50 


3.67 


22 


0.52 


0.26 


0.65 



V254. Assuming that the spectrum originates from the sec- 
ondary leads once more to amplitude ratios AV c /AV i, s which 
for all mass ratios are much smaller than unity. Upon adopting 
that the source of the spectrum is the primary we get the results 
displayed in Table [5J which suggests that V254 is similar to 
V214. 



Table 5. Derived parameters for V254 



1 


a 


i 


mi 


WJ2 


AVJAV obs 




[Re] 


n 


[M e ] 


[Me] 




0.2 


2.77 


72 


1.61 


0.32 


1.7 


0.3 


2.98 


43 


1.84 


0.55 


1.0* 


0.5 


3.29 


24 


2.15 


1.07 


0.6 



NV334. Because of the low brightness of the system 
(V = 18.83 mag) we were able to obtain only a very 
rough estimate of the velocity amplitude. We therefore modifed 
our analysis: we estimated the inclination based on the light 
curve, and compared the calculated velocity amplitude K c to 
the observed one. Assuming that the spectrum originates in the 
secondary leads to K c being much larger than K a b s , regardless 
of the value of q. The same analysis applied to the configuration 
with the primary generating the spectrum produces results 
shown in Tabled which seem to indicate that NV334 is similar 
to V251. However, the similarity is superficial: the effective 
temperatures of the two systems differ by ~ 4000 K, causing 
the cooler and dimmer NV334 to occupy an entirely different 
position in the color-magnitude diagram of cj Cen (see Fig.Q}. 



Table 6. Derived parameters for NV334 



1 


a 


i 


mi 


m.2 


K c /K obs 




[Re] 


n 


[Me] 


We] 




0.1 


1.22 


76 


0.33 


0.03 


0.4 


0.3 


1.45 


65 


0.47 


0.14 


1.2* 


0.5 


1.60 


63 


0.55 


0.27 


1.8 



NV360. Our standard analysis fails for this system (iterations of 
the inclination i based on the velocity curve are either converg- 
ing extremely slowly or entirely diverging), so that we again 
have to apply a modified version based on the light curve. As for 
the case of NV334, assuming that the spectrum originates in the 
secondary leads to K c being much larger than K t, s , regardless of 
the value of q. Systems with the primary being responsible for 
the spectrum fare much better (see Table [7), and we conclude 
that NV360, unlike NV334, is truly similar to NV251. 



Table 7. Derived parameters for NV360 



1 


a 


i 


mi 


ffl2 


K c IK obs 




[Re] 


[°] 


[Me] 


[M s ] 




0.5 


2.52 


26 


0.36 


0.18 


0.7 


0.8 


2.79 


25 


0.41 


0.32 


1.0* 


0.9 


2.86 


25 


0.42 


0.37 


1.1 



NV400. When the secondary is assumed to generate the spec- 
trum, the calculated light-curve amplitude is for all values of q 
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much higher than the observed one. Assuming that the spectrum 
originates in the primary we get the results displayed in Table [8] 
Apparently, N V400 is yet another system similar to V25 1 . 



Table 8. Derived parameters for NV400 



q 


a 




mi 


Wi2 


AV c /AV obs 




[R e ] 


[°] 


[M ] 


[M G ] 




0.3 


2.65 


28 


0.48 


0.14 


1.3 


0.5 


2.93 


13 


0.56 


0.28 


0.8* 


0.8 


3.24 


9 


0.63 


0.50 


0.2 



V240. Assuming that the spectrum originates from the primary 
we obtain the results shown in Table [9] computed light-curve 
amplitudes are far too low (if the secondary were the source of 
the spectrum they would be even lower). Moreover, the total 
mass of the system is unacceptably high (as we argued in the 
case of V214, it should not exceed 1.84 M ). In other words, the 
color and apparent magnitude of V240 are incompatible with its 
observed velocity amplitude, suggesting that this system does 
not belong to to Cen. Another possibility is that the light of 
V240 is dominated by a tertiary component, and the recorded 
spectra are unrelated to the photometric binary. 



Table 9. Derived parameters for V240 




-0.5 0.5 

phase 




Fig. 4. Light and velocity curves of blue stragglers - continua- 
tion of Fig. [3] A single spot 2% hotter than the photosphere was 
placed on the primary of V20 in order to qualitatively account 
for the asymmetry of the light curve (just for illustrative pur- 
poses; it was not included while analyzing the data). Parameters 
of the fits for NV400 and V20 are given in Tables[8]and[n]in the 
rows indicated with an asterisk. No fit was possible for V240, 
because the velocity ampliutude turned out tobe smaller than the 
observational errors. The light curve, however, is coherent and 
remained stable between 1999 and 2003. 



<? 


a 


i 


nti 


m.2 


AV c /AV Bbs 




[Re] 


[°] 


[M s ] 


[Me] 




0.07 


2.83 


16 


2.59 


0.18 


0.1 


0.1 


2.96 


11 


2.88 


0.29 


0.05 






-0.5 0.5 

phase 




Fig. 3. Light and velocity curves of blue stragglers. Parameters 
of the plotted fits are given in Tables [3] |4] [5] and [7] in the rows 
indicated with an asterisk. 




0.5 0.5 

phase 




0.5 

phase 



Fig. 5. Light and velocity curves of subdwarfs. Parameters of the 
plotted fits are given in Tables l6landfT2lin the rows indicated with 
an asterisk. 



3.3.2. NGC 6397 

V17. The visible component must be the secondary (an alterna- 
tive assumption enforces inclinations at which, regardless of the 
value of q, the ratio AV c /AV h s is much lower than unity). Our 
standard analysis produced the results shown in Table [TOl 

The fact that V17 is situated at the turnoff of the main se- 
quence prompted us to consider an exceptional case of a de- 
tached binary with mi e qual to the turnoff m ass, which in NGC 
6397 is close to 0.8 M Q dKaluznv et al.f2 008). Keeping a conser- 
vative approach we set ma = 0.75 M Q , and, as expected based 
on Table [TOl we found that the amplitudes of all calculated light 
curves were too low. Our results indicate that the primary, al- 
though dimmer than the secondary, is rather massive, with a 
mass possibly exceeding 1 M . 
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-0.5 0.5 -0.5 0.5 

phase phase 



Fig. 6. Light and velocity curves of systems at the main se- 
quence turnoff. Small ticks in light-curve plots are separated by 
0.01 mag. Parameters of the plotted fits are given in Tables [TUl 
and[T3]in the rows indicated with an asterisk. 

Table 10. Derived parameters for V17 



9 


a 




m 1 


Wl2 


AV c /AV obs 




[Re] 


n 


[M ] 


[M G ] 




0.3 


6.00 


14 


1.99 


0.60 


0.7 


0.4 


5.56 


16 


1.47 


0.59 


0.9 


0.5 


5.25 


19 


1.15 


0.58 


1.1* 


0.6 


5.02 


21 


0.95 


0.57 


1.5 



V20. This is the only system in our sample for which it was pos- 
sible to measure the velocities of both components. The mea- 
sureme nts were performed usin g the TODCOR method devel- 
oped bv lZucker & Mazehl(ll994l) . which also calculates the ratio 
of component fluxes a = fi/fi- Based on all available spectra of 
V20 (there were 9 of them), we obtained a = 2.47 + 0.02 which, 
combined with M" , = 2.94 mag, yields M° b n = 3.31 mag and 
M° boa = 4.29 mag" 

Since the mass ratio is small (q « 0.2), the secondary must 
be rather oversized for its mass. We assume that it fills its Roche 
lobe, and we apply the following nonstandard procedure: i) set i; 

ii) iterate a and q to fit the velocity curves of both components); 

iii) adjust the size of the primary to get M boll = 3.31 mag; adjust 
the temperature of the secondary T% to get M bgl2 = 4.29 mag 
(note that this is the minimum possible temperature). The best 
fit among those displayed in Table QT| indicates that V20 is 
another system similar to V25 1 . 



Table 11. Derived parameters for V20 



i 


a 


1 


mi 


m.2 


AVJAV obs 




[Re] 


[°] 


We] 


[Me] 




30 


5.28 


0.20 


2.23 


0.45 


0.4 


40 


4.09 


0.20 


1.04 


0.20 


0.8 


45 


3.72 


0.20 


0.78 


0.15 


1.0* 


50 


3.43 


0.20 


0.61 


0.12 


1.2 



V27. Upon applying the standard analysis we found that, re- 
gardless of whether the primary or the secondary produces the 
spectrum, for a given value of q either there is no inclination 
for which the calculated amplitude matches the observed one or 
the masses of both stars are unacceptably low {m\ < 0.05 M Q ; 
ni2 < 0.04 M Q ). One can make them higher by assuming that the 



system is detached instead of semi-detached; however, in such a 
case it is impossible to find firm lower mass limits. 



Table 12. Derived parameters for V27 



q 


a 




nt\ 


m.2 


T 2 


AVJAV obs 




[Re] 


[°] 


[Me] 


[Me] 






0.2 


2.00 


55 


0.29 


0.06 


6500 


1.0* 


0.2 


2.50 


40 


0.57 


0.11 


6000 


0.9 


0.2 


3.00 


33 


0.98 


0.20 


6000 


0.3 


0.3 


1.80 


43 


0.20 


0.06 


6000 


0.9 


0.3 


2.00 


36 


0.27 


0.08 


6000 


0.7 


0.3 


2.50 


28 


0.52 


0.17 


6000 


0.2 



To obtain a sample of detached solutions we modified our 
standard procedure in the following way: i) set q and a; ii) adjust 
i so that the calculated velocity curve agrees with the observed 
one; iii) adjust the size of the primary to get M\, \\ = 6.0 mag; 
iv) adjust the size (and, if necessary, the temperature) of the sec- 
ondary to get Mboii = 6.3 mag (the latter two values are of course 
arbitrary, the only constraints being that the secondary is not 
much dimmer than the primary, and the combined bolometric 
magnitude of the system is 5.38 mag), v) check if the calculated 
light curve agrees with the observed one. 

If one assumes that the spectrum is generated by the sec- 
ondary, then the inclination consistent with the velocity curve is 
much too low to reproduce the light curve (except for systems 
with q > 0.75, but in those cases the masses of the components 
are smaller than 0.05 M ). The only reasonable solutions we 
were able to find for the detached configuration are based on the 
assumption that the primary generates the observed spectrum, 
while the secondary is mainly responsible for the observed 
ellipticity effect. The results shown in Table fT2l seem to favor a 
configuration with small and rather discrepant masses. 

V36. The visible component must be the secondary; the alter- 
native assumption enforces inclinations such that the calculated 
amplitude of the light curve is much lower than the observed 
one for all values of q. However, assuming that the secondary 
fills its Roche lobe is untenable, as our standard analysis indi- 
cates that in such a case the calculated amplitudes of the light 
curve are too high. Thus, we are forced to adopt a detached con- 
figuration. Since V36 is located at the turnoff (only 0.3 mag be- 
low V17; see Fig. we assume that mi - 0.75 M and apply 
the following four-step procedure: i) set q\ ii) adjust a to get 
m.2 = 0.75 M ; iii) adjust i so that the calculated velocity curve 
agrees with the observed one; iv) adjust the size of the secondary 
to get Mbou = 3.65 mag. The results shown in Table [T3l indicate 
that V36 is similar to V17, with the primary's mass possibly even 
higher than 1.5 M Q . 



Table 13. Derived parameters for V36 



1 


a 


i 


m.2 


m x 


AV c /AV obs 




[Re] 


[°] 


[M~e] 


We] 




0.2 


7.40 


20 


0.75 


3.11 


0.9 


0.3 


6.62 


24 


0.75 


2.49 


1.0* 


0.4 


6.17 


28 


0.75 


1.87 


1.3 


0.5 


5.86 


32 


0.75 


1.50 


2.1 
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4. Discussion and conclusions 

10 out of 1 1 objects in our sample exhibit radial-velocity varia- 
tions indicating their binary nature (the only exception is V240 
in u> Cen). In that sense, the sample was well chosen. Regarding 
the principal aim of the present survey we have been less suc- 
cessful - no clear-cut evidence for high-mass degenerate compo- 
nents was found. However, while 8 systems proved to be more 
or less ordinary binaries, the remainig two (V17 and V36 at 
the turnoff of NGC 6397) cearly deserve further scrutiny. First, 
they are the only systems in which the brighter component is 
the less massive secondary. Second, the masses of their dim pri- 
maries may be significantly larger than 1 M G (in V36, even larger 
than 2 M Q if the secondary is indeed a turnoff-mass star). Third, 
they are weak X-ray s ources. Based on Chandra observations, 
iBogdanov et al.1 d2010h classify them as active binaries (AB), 
i.e. systems composed of main-sequence or subgiant stars whose 
weak X-ray emission originates from magnetic activity. Our re- 
sults rule out the possibility that they are composed of pristine 
stars which have not undergone any mass transfer episode, as 
in such a case the more massive compo nent would also have to 
be the brighter one. This is not to say iBogdanov et al ] (1201 Oh 
are wrong: short-period degenerate binaries in quiescence may 
easily "masquerade" as AB systems, since they are most likely 
synchronized and their visible components spin fast enough for 
the stellar dynamo to be highly efficient. 

According to the most optimistic interpretation of our re- 
sults, V17 and V36 may contain a neutron star (the dim pri- 
mary of V36 may even be a black hole). Such systems can in- 
deed be expected, as population synthesis calculations indicate 
that the formation of a degenerate binary in which a neutron 
star or a black hole is accompanied by a main-sequence star 
or a subgiant is nothing unusual in globular clusters. Although 
the rate at which they form is not high (< 2.5 systems/Gyr), 
they should be transient all the time, and therefo re more likely 
seen as a qLMXBs rather than as bright LMXBs dlvanova et alJ 
[2QQ8h . Moreover, NGC 6397 bel ongs to core-collapse clusters 
for which the formation rates of llvanova et al.l (120081) are just 
lower limits. Another (admittedly, not too strong) support for the 
optimism comes fro m the detection of a possible qLMXB near 
the turnoff of M30 dLugger et al.l [2006). Main-sequence secon- 
daries similar to those in V17 and V36 are also found in galactic 
X-ray binaries whose primary components are black hole candi- 
dates; e.g. V1033 Sco, GRS 1739-278 or V821 Ara (Ziolkowski, 
priv. comm). 

The fact that all our blue stragglers except V240 are binary 
supports the currently leading hypothesis concerning the nature 
of these objects, namely that they result rather from an extensive 
mass exch ange between the c omponent stars than from stellar 
collisions dKnigge et al.ll2009h . In all cases the brighter compo- 
nent is the primary which according to the mass-transfer sce- 
nario must have acquired significant amounts of hydrogen-rich 
material from the envelope of the originally more massive sec- 
ondary. In fact, our results seem to indicate that blue stragglers 
are Algol-like systems in which the original mass ratio has been 
reversed, causing the mass transfer to effectively stop. Our strag- 
lers may be similar to the well-studied star V228 in 47 Tuc 
Kaluznv et aT1l2007l) . however we cannot exclude the possibil- 
ity that in some of them the primary is already at the beginning 
of the subgiant branch and has become large enough to approach 
its Roche lobe. 

V214 and V254 in to Cen have primaries with masses likely 
exceeding 1 M (maybe even 1.5 M ), and their total masses 
approach 2 M . Systems that massive must have originated 



from binaries subject to conservative or nearly conservative mass 
transfer. They also have large amplitudes of light variations 
(0.12, mag and 0.10 mag, respectively), in our sample rivaled 
only by that of NV334 (~0.3 mag); however a significant contri- 
bution to this effect must originate from relatively large inclina- 
tions found for all three systems. The remaining blue stragglers 
are significantly less massive - their total masses are lower than 
1 M . As far as we know, our results are the first to indicate such 
a broad (and possibly bimodal) mass distribution. If this effect 
is real, it probably reflects differences in the efficiency of mass 
transfer and/or mass loss. We note in passing that some mass 
may still be flowing between the components of V20, as this 
system exhibits an asymmetric and variable light curve, and is 
a weak X-ray source. Another possibility to account for the ap- 
parent large width of the mass distribution is related to the pecu- 
liarities of the chemical composition of a> Cen, in which strongly 
He-enriched subpop ulations are suggested to exist, w ith Y reach- 
ing up to ~ 4 (e.g , lNorrisl2004HPiotto et aT1l2005h but see also 
ICatelan et al.l 120 1 Ol for a recent review and some caveats). Our 
fits are based on the standard value Y = 0.24, and for a given 
B - V color they can overestimate t he temperature of e nriched 
stars by up to a few hundred Kelvin dDotter et al]l2008h . If this 
is what happened, then our standard analysis caused the He-rich 
systems to artificially shrink both in size and mass. Such hypoth- 
esis would also explain why stragglers with signficantly different 
"apparent" masses are able to form a relatively tight group on the 
H-R diagram. 

The subdwarfs NV334 and V27 are very similar to each 
other in that they have large mass ratios and low-mass primaries 
with practically the same effective temperatures. Our estimate of 
the primary's mass in NV334 is consistent with the canonical hot 
subdwarf mass of 0.47 M , while the estimated mass of the sec- 
ondary falls near the peak of the mass dis tribution of the un seen 
companions to field hot subdwarf stars ( Geier et ail 20091). O n 
the other hand, and as pointed out by Moni Bidin et al. 1 120081) . 
the close binary fraction among subdwarf stars appears to be 
much lower in GCs than in the field, which may point to different 
formation mechanisms for many of the GC subdwarfs. In V27 
the estimated masses are significantly smaller (see Table [T2l ; 
note however that they are also less reliable due to problems dis- 
cussed in Section [3.3.21 and that the solution with m\ = 0.57 M Q 
and 1112 = 0.11 M Q is also acceptable. Taken at face value, the 
starred parameters in Table [12] would indicate that V27 is very 
similar to HS 2231+2241, an HW Vir-type system with a com- 
panion at or below the very low mass limit for M-dwarfs. The 
~ 0.26 M primar y of HS 2231+2241 i s a non-helium burn- 
ing, post-RGB star (Os tensen et alJ[2008h . Another example of 
such a type of sta r is the 0.24 M component of HD 188112 
dHeber et alj|2003l) . The primaries of NV334 and V27 are much 
cooler (T ~ 7500 K compared to T = 28400 K and T = 20500 
K, respectively, for HS 223 1 +224 1 and HD 1 88 1 1 2), but this dif- 
ference may result from their more advanced evolutionary stage. 

Our sample is too small for a meaningful statistical analysis; 
nevertheless the above-discussed results indicate that further re- 
search on ellipsoidal binaries in globular clusters is a wortwhile 
task, and suggest its most promising direction: it seems that fu- 
ture sampling of the photometric variables should be biased to- 
wards objects located on or to the right of the main sequence, 
including those identified as weak X-ray sources. Obviously, for 
purely economical reasons (exposure times!) the targets sould 
not fall much below the turnoff point. Such an approach is te- 
dious, but at the same time it is the only one on which a fair 
census of degenerate binaries in globular clusters can be based. 
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